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There are different mechanisms proposed to be responsible for the temperature dependence of 
the interlayer exchange coupling (IEC), namely a smearing out of the spacer or interface Fermi 
surface and excitations and interactions of spin waves. We propose a possibility to separate both 
effects by calculating the excitation spectrum of an extended Heisenberg model and connecting the 
results with ferromagnetic resonance (FMR) experiments. To solve the Heisenberg model we use 
an approximation that was shown to yield excellent results. In this paper the main idea of this 
procedure is explained and a detailed investigation of the spin wave contribution to the temperature 
dependence of FMR resonance frequencies and fields is carried out. 



I. INTRODUCTION 

The discovery of the interlayer exchange coupling 
(IEC) has triggered a renewed interest in the physics of 
ultrathin metallic films in the last decades. The term 
IEC denotes the coupling of magnetic layers separated 
by non magnetic spacer layers that occurs in quite a 
number of metallic and semiconducting systemsi* The 
general mechanism of the coupling and a lot of addi- 
tional features are well understood today. The coupling 
is caused by spin dependent reflection of spacer electrons 
at the magnet/non-magnet interfaces leading to a spin 
dependent interference and therewith to a spin depen- 
dent renormalization of the density of states (DOS) and 
the free energy within the spaceri The modification of 
the spacer DOS and energy depends sensitively on the 
relative orientation of the magnetic layers favoring a spe- 
cial configuration. This results in the interlayer coupling. 
The coupling may be ferro- or antiferromagnetic and os- 
cillates with respect to the spacer thickness. The periods 
of the oscillation depend on the spacer Fermi surface, 
while the amplitudes as well as the phases may be addi- 
tionally influenced by interface roughness, interdiffusion, 
disorder, and related effects. 

However, up to now it is not clear which kind of pro- 
cess dominates the temperature dependence of the cou- 
pling. There are several candidates which are evaluated 
and compared in Ref . 

(i) One source of temperature dependence is con- 
nected with the described coupling mechanism and 
is caused by the smearing out of the Fermi-surface 
of the spacerA^ Furthermore temperature influ- 
ences the reflection coefficients at the magnet/non- 
magnet interfaces £ Both processes influence the 
energy of the spacer layer and will be called "spacer 
effect" in the following. 

(ii) Another source of temperature dependence are 
magnetic excitations, altering the properties of the 
magnetic layers^ ("magnetic layer effect"). Ther- 
mal magnetic disorder may drastically reduce the 
energy difference between parallel and antiparallel 
alignment of the magnetic layers and therewith the 



interlayer coupling. 

To find the dominating mechanism elaborate ferromag- 
netic resonance (FMR) measurements were performed at 
a prototype trilayer^ namely a Ni7/Cu x /Co2/Cu(001) 
system where the subscripts denote the number of atomic 
layers of each film. The data are analyzed using a clas- 
sical continuum model, where the classical equation of 
motion for the magnetization 10 or an expansion of the 
free energjiii is considered. This analysis gives the tem- 
perature dependence of an effective interlayer coupling 
J/(T) which is discussed in detail in Ref. [12 and recon- 
sidered in Ref. [jj This quantity is proportional to the 
energy difference between the state where the magneti- 
zation of both magnetic films are parallel to each other 
and the state where the magnetization of the lower mag- 
netic film are aligned antiparallel to the magnetization of 
the upper one: 

j T (T) ~ F^(T) - F n (T). (1) 

Note that the magnetic moments are supposed to be par- 
allel within the magnetic films. The most important fea- 
ture of JiiT) found experimentally is its effective T 3 / 2 - 
like decrease with temperature. Naturally both tempera- 
ture dependencies, i.e. the temperature dependence due 
to the Fermi surface softening and those due to mag- 
netic excitations, are included in Jj(T). Thus, since both 
mechanisms (i) and (ii) contribute, the dominating mech- 
anism cannot be deduced directly from Jj(T)miH 
On the other hand a procedure sensitive to only one of 
the mechanisms (i) or (ii) would be highly desirable. It 
would allow to separate the mechanisms and to deduce 
the source of temperature dependence of the IEC directly 
from the measurements. It is the intention of this paper 
to show that such a procedure exists. The new approach 
is based on an alternative analysis of FMR measurements 
using the microscopic Heisenberg model. 
To illustrate the idea let us consider a modeling of the 
considered trilayer systems: 

(ij)a (iajp) 
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(2) 
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where a runs over the two different layers while i and 
j denote a lattice site within a layer. For convenience 
the sums in the first lines are taken over nearest neigh- 
bors only. Thus the iniralayer coupling parameters J a 
as well as the interlayer coupling parameter J/ are ef- 
fective parameters containing also the contributions of 
non-nearest-neighbor interactions. The model describes 
two magnetic layers coupled by the interlayer coupling 
Jj, which may be positive or negative. The intralayer 
coupling constants J a are chosen to be positive to de- 
scribe ferromagnetic layers. The model further accounts 
for anisotropies due to spin orbit (K201) and dipolar (go a ) 
coupling as well as for an external field gj Q /isBo applied 
in arbitrary direction. The parameter gj a denotes the 
Lande factor of the layer a. The anisotropy due to spin- 
orbit coupling K2a, which includes magnetocrystalline 
as well as magnetoelastic contributions^ will often be 
termed "lattice anisotropy" in the following. The film 
plane is the xy-plane and tetragonal symmetry is as- 
sumed. Such a model is often used to describe mag- 
netic properties of transition metal films^* especially in 
the presence of anisotropies. Note that the spacer of the 
trilayer system is not considered explicitely but its ef- 
fects on the magnetic layers are described by the model 
parameter J/. The latter is proportional to the energy 
difference between the states of parallel and antiparallcl 
aligned magnetizations in both layers, M || M& and 
M a fj, Mb, if saturated magnetism is assumed. This ap- 
plies for T = or at finite temperatures for cases where 
the spin wave excitations of the model are artificially 
frozen. For the latter case a temperature dependence of 
J/ can be caused exclusively by excitations within the 
spacer. Hence, Ji{T) is proportional to the difference 
between the energy contributions of the spacer in the par- 
allel and antiparallel state: 

Ht)~f;v{t)~f;i{t). (3) 

Now we want to connect the model @ with the men- 
tioned FMR experiments. The probed quantity of an 
FMR experiment is the so called "resonance frequency" 
that is identical to the spin wave mode with infinite wave- 
length uio& This mode is also called uniform (q = 0) 
mode. In a coupled bilayer system there are two such ex- 
citations, termed optical and acoustical mode (wo ,^0a)- 
This can easily be understood if a semiclassical spin 
model is considered. The uniform mode of a mono- 
layer corresponds to a spin wave where all spins rotate 
in phase. The acoustical mode of the coupled bilayer 
system corresponds to the state where the spins in both 
layers rotate in phase. If the optical mode is excited, the 
spins in the upper layer rotate in anti-phase compared to 



the spins of the lower layer. Depending on the coupling 
sign, the optical mode can be found at higher (ferromag- 
netic interlayer coupling) or at lower frequencies (anti- 
ferromagnetic coupling). Note that the FMR technique 
has equal access to ferromagnetic and antiferromagnetic 
coupling. In this picture the distance between the modes 
Aco>o is a measure for the interlayer coupling strength if 
equivalent magnetic layers are assumed. 
In practice, during an FMR experiment the probe fre- 
quency Vhf is fixed in the microwave region, while the 
magnitude of the external field is varied until the reso- 
nance frequency equals the probe frequency. 

2irv hf = wo( ,o) ( B o) (4) 

This relation is fulfilled at special magnitudes of the ex- 
ternal field called "resonance fields" (-E>o(a,o))- The val- 
ues of the resonance fields as functions of the direction 
of the applied external fields (-B Ot .( O)O )(0-B o )) ^ s tne mos t 
important experimental output^ For an appropriately 
chosen probe frequency there are also two resonance fields 
in a coupled bilayer system belonging to the optical 
and acoustical mode. Again, the distance between the 
resonance fields (ABo r ) is influenced by the interlayer 
coupling strength J/ and the anisotropies in both layers. 
An FMR experiment can be analyzed using the model (0 
by calculating the spin wave modes Wo(a,f>) an d the res- 
onance fields -Bo(a,o) and fitting the model parameters 
until the calculated values correspond to the measured 
ones. It is possible to obtain the interlayer coupling as 
well as the anisotropy strengths from the angular depen- 
dence of the resonance fields. 

Let us now consider the temperature dependence of the 
resonance frequencies and fields. There are two principal 
sources of temperature dependence: 

1. If the parameters of the model J5J are fitted to real 
systems as proposed above, they may be tempera- 
ture dependent itself. This applies especially to the 
interlayer coupling parameter Jj which is not a fun- 
damental microscopic coupling constant but rather 
given in Eq. @ and which may be significantly 
influenced by temperature As seen in Eq. @ 
the temperature dependence of Ji(T) describes the 
temperature dependence of the spacer energies. 

2. Additionally, even if the model parameters are 
fixed, there is an intrinsic temperature dependence 
of u>o(a,o) and A_B ( a , o ) m the model J2J) due to inter- 
actions of spin waves excited within the magnetic 
layers. 

It is obvious that the temperature dependence of the 
model parameter Ji(T) can be identified to the "spacer 
effect" (i) while the intrinsic temperature dependence due 
to spin wave interactions is equivalent to the "magnetic 
layer effect" (ii). 

In this situation the microscopic model J2J ca n be used 
to study the mechanisms (i) and (ii) separately. On the 
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one hand one can fit the model parameters to experimen- 
tal FMR data which gives the function Jj (T) . The latter 
exclusively describes the temperature dependence of the 
IEC due to the "spacer effect" (i). 

On the other hand the model parameters and especially 
Ji can be fixed. Then one can investigate the tempera- 
ture dependence of the resonance frequencies due to spin 
wave interactions. Thus the temperature dependence of 
the FMR signal caused by the magnetic layer effect (ii) 
alone can be analyzed. 

The motivation of our paper is twofold. First we want 
to show that an evaluation of FMR data is in principle 
possible by use of the microscopic model However, 
our paper is mainly dedicated to the investigation of the 
second point, i.e. of the temperature dependence of the 
FMR signal caused by the interaction of spin waves. 
The paper is organized in the following way: 
Using our microscopic model we will discuss two quanti- 
ties: the resonance frequencies u>o(a,o) an d the resonance 
fields So r ( a>0 ). We will consider a monolayer (simulat- 
ing a single film) as well as a system of two coupled 
layers (which simulates prototype IEC-trilayer systems). 
We compare the results of the classical model with those 
of the microscopic model at zero temperature before we 
discuss the intrinsic temperature dependence of the res- 
onance frequencies and fields in detail. With this proce- 
dure we want to answer the following central questions: 

• Is there a relevant intrinsic temperature depen- 
dence of the uniform spin wave modes in the model 
iPJ and how is it characterized? 

• What are the differences between the results of our 
microscopic approach and of the macroscopic clas- 
sical model for zero temperature, in other words, 
can the microscopic model give new insight in FMR 
experiments at a fixed temperature, too? 

It turns out that the latter is indeed the case. For a 
nickel films grown on a Cu(OOl) surface the absence of a 
resonance field in the hard direction was reported— This 
observation follows directly for larger out-of-plane lattice 
anisotropies in the microscopic model, while it can be 
hardly explained within the classical model. Furthermore 
we will find a significant temperature dependence of the 
uniform modes due to spin wave interactions. A descrip- 
tion of this temperature dependence in terms of effective 
model parameters (K2(T) and J/(T)) is possible approx- 
imately only and justified only for lower temperatures 
T <T C . Before turning to the presentation of our results 
we will shortly discuss the approximations used to solve 
the model iffil. The intended use poses clear requirements 
to the theory. Its numerical effort must be moderate to 
keep the fitting procedure feasible, while it has to yield 
accurate results for the spin waves for arbritrary direc- 
tions and magnitudes of the external field. Such a theory 
is proposed and evaluated in Ref. and applied to our 
model (J2J) in the next section. As already mentioned this 
contribution uses results and insights from earlier works, 



namely from Ref. and Ref. E The former (Ref. 
0) discusses and compares the different mechanisms that 
may cause a temperature dependence of the interlayer ex- 
change coupling on a very general level therewith clearly 
formulating the problem we want to consider in this pa- 
per. The latter (Ref. HI) is mere technically proposing a 
theory to solve our model @- 



II. THEORY 

In this section we want to motivate the theoretical ap- 
proach used to solve the extended Heisenberg model J2J- 
Due to the tetragonal symmetry the directions in the film 
plane are equivalent. This allows us to restrict our calcu- 
lations to a plane perpendicular to the film, the xz-plane. 
The external field and the magnetization are confined to 
this plane in the following: 

Bo = {Box, 0, B 0z ) 
(S) = ((S x ),0,(S z )). 

To proceed we split the Hamiltonian J3J) into two parts. 
The first one, 

Hi = S J a Si a Sj a S JjSi a Sjfj 

(ij)a {iajfi) 



(5) 



contains all terms except for the dipolar coupling. In a 
previous paper— we developed a theory for the Hamil- 
tonian JSJ for the special case of a monolayer. It 
was shown by comparison with Quantum Monte Carlo 
(QMC) results^ that our approximation is in fact very 
accurate for arbritrary directions of the external field 
and can therefore be used here. Our proposal combines 
an RPA decoupling procedure 20 applied to the nonlocal 
Heisenberg exchange terms and a generalized Anderson 
Callen approximationSLSii for the local anisotropy term. 
Details are given in Ref. The theory is based on a 
transformation into a local coordinate system £ — » £' 
where the new z'-axis is parallel to the magnetization. It 
is easy to generalize this theory to a multilayer system if 
the RPA decoupling is also used for the interlayer cou- 
pling term. A remarkable result of this theory is that the 
anisotropy and the interlayer coupling can be replaced by 
internal fields Bi Q in the single magnon Green function 
matrix of the transformed system: 



(6) 



The effective field is a sum of the external field, a term 
due to the anisotropy and a further term resulting from 
the interlayer coupling: 



Blc — JJa^lflBo + Bk 2 a + B 



J I a • 



(7) 
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The interlayer field is given by 

B Jia = 2J 7P '(S 7 ), 



(8) 



where p' denotes the number of nearest neighbors in the 
adjacent layer 7. For the anisotropy field one finds 



B 



K2 Ct X 



-K 2a (S ax )sm 2 6 a C{ 



B K2az = +2K 2a (S az ) i 1 - - sin 2 9 a j C[ 
with 

C'i = l-^{S{S + l)-{Sl z ,)), (9) 

where 9 a is the angle between the z-axis and the magne- 
tization in layer a. Note that the anisotropy field B/f 2 is 
not an analytic function of the magnetization (S z >) only, 
but of (5 2 /)-terms as well. The single- magnon Green 
function 10 is now a functional of the effective fields: 



{ (2S bz >) 



via the matrix M r 



(2J a (S az '){p - 7 q ||) - B la ) 
(2Ji(5 /9 » / >7q± cos(0 a -0 6 )) 

(Mfy. 



(10) 



Generally, the magnetization angles in both layers, 9 a 
and 9b, may be different. Note that the coordinate system 
is rotated by 9 a in layer A and by 9b in layer B. Thus the 
z'-axes in both layers are not parallel to each other. The 
respective angles are determined by the effective fields 



tantL = 



B n 



B n 



(11) 



Eq. I|1U|> is valid as long as the deviations between the 
magnetization angles of different layers are small. 
The coordination number within a layer is denoted by p 
and the structural factors 7 q || and 7 q j_ are a consequence 
of the two dimensional Fourier transformation and de- 
pend on the lattice structure. For the fcc(100) geometry 
considered here they are given by 

1 1 

7 q || = 4 cos -aq x cos -aq y 



7q± = [2cos^aq y + e~ l i aq * 



where a is the lattice constant. 

From the single-magnon Green function (|1U[) one can eas- 
ily obtain the expectation values (S z i) and (Si,) using the 
standard text book procedure 2 ' which can be generalized 
to film geometries straightforwardly, 2 - Additionally the 
spin wave frequencies hvq = are directly given by 
the eigenvalues i? q of the matrix M q . Eqs. 0) - 
establish a self consistent set of equations and thus the 



problem is solved as far as the Hamiltonian (|3J is con- 
cerned. 

Now the dipolar term 
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S?aS^ a — (Si a Yij ) (Sj a Vij ) J (12) 



has to be treated. The interaction is nonlocal and very 
similar to the Heisenberg exchange terms. Thus an RPA 
decoupling would be certainly a good choice. However it 
was shown in Ref. |2{| that the more feasible mean field 
decoupling also yields acceptable results. A further term 
has to be added to the effective field due to the dipolar 
coupling: 



where the dipolar field is given by 2 ^ 

Bdip a x — ffOa (Sxa ) ^ ^ 3"" 



(13) 



'Oj 



-Bdipaz — 2gQ a (S za ) ^ ' 

i T 0j 



(14) 



The sum runs over all lattice sites of a layer. There are 
no further changes in the relations Q- Ijllfl due to the 
dipolar coupling. Thus the effective field B Q defined in 
Eqs. (|13fl and (|14(l simply has to replace the field Bi Q 
from Eq. J7J if one wants to account for dipolar cou- 
pling. 

As already mentioned before, this theory is capable of 
calculating the whole spin wave excitation spectrum and 
especially the uniform modes u)q and tooa- Therewith 
the resonance fields B 0r ^ a o ^ can be obtained for any mi- 
crowave probe frequency v^f. This holds for arbritrary 
external field directions and for any temperature. Thus 
one can calculate the central functions B 0r ( a :O )(9b ) and 
Bo r (a,o)(T) as well as the dependence of the resonance 
frequency on the external field uio^ a>0 )(B). The latter 
is often referred to as "dispersion relation" in connec- 
tion with FMR experiments and gives the excitation gap 
-E-gap = fr-nrin(cL>( 0a ), w (0o))- The results will be presented 
and discussed in the next section. 



III. THE SINGLE FILM 

In the following sections the anisotropy parameters and 
the dipolar coupling will be given in units of [/iskG] with 
1 /iskG = 5.79 /ieV. As can be seen in the Hamilto- 
nian J2J the parameters Jj, J ai K 2 and go denote energies 
per lattice site, i.e. energy densities rather than energies. 
The external fields are given in [kG] while the resonance 
frequencies v§ — luq/2tt are given in [GHz]. The spin 
quantum number is set to unity (S=l) and the Lande 
factor is fixed to the bulk nickel value gj a = 2.21. 
Though we are mainly interested in the coupled bilayer 
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FIG. 1: The resonance frequency as a function of the ex- 
ternal field at T — for different angles between the easy 
direction and the external field A8b (symbols). Left panel: 
positive lattice anisotropy Ki = 10 ^kG, right panel: neg- 
ative lattice anisotropy K2 = — lO^kC, inset: dipolar cou- 
pling go — 2.1/iskG. Here and in the following pictures the 
spin quantum number is set to unity (S=l). The results of 
the classical model are also shown (solid lines). 




B (kG) 

FIG. 2: The dispersion relations vo(Bo) at four different re- 
duced temperatures t = T/T c . The external field is applied 
in the easy (upper curves) and in the hard direction (lower 
curves). The lattice anisotropy is K2 = 15 /iskG and the 
dipolar coupling is zero. The thin solid line gives the reso- 
nance frequencies for paramagnetism. Circles: t = 0, K2 = 
8.2 ^ s kG 



system one has to understand the single film system be- 
fore a meaningful interpretation of more complex systems 
is possible. Moreover, though the experimental output is 
the resonance field it is very instructive to study the "dis- 
persion relation" . The dispersion relations for zero tem- 
perature are shown in Fig. ^ The spin wave frequency Vq 
is displayed as a function of the external field B for pos- 
itive and negative lattice anisotropics K 2 = ±10/iskG 
(left and right panel, go = 0) as well as for dipolar cou- 
pling g = 2.1/iskG (inset, K 2 = 0). The angle A9b 
denotes the angle between the external field and the easy 
direction. The latter is perpendicular to the film for a 
positive lattice anisotropy and parallel to the film plane 
for negative lattice anisotropics or for dipolar coupling. 
For comparison the results of the classical model as e.g. 
given in Ref. [3 are also shown (thin solid lines). The 
anisotropy parameter K2 used for the classical calcula- 
tion is scaled as proposed in Ref. 26 (K 2 = S(S — l/2)K 2 ) 
and accounts for both, lattice and dipolar anisotropies. 
Qualitatively the results are similar for all cases. If the 
external field is applied parallel to the easy direction 
(AS b = 0°) the function vq(Bq) increases linearly. The 
intersection with the ordinate gives the minimal excita- 
tion gap due to the anisotropy. This gap is larger for 
positive than for negative lattice anisotropies and rises 
with the norm of K 2 as can be seen by a comparison of 
the left panel of Fig. Q (K 2 — lO/i^kG) and the solid 
line in Fig. |2 (if 2 — 15/iskG). There are more features 
in the dispersion relation if the external field is applied in 
the hard direction (A0g o = 90°). The magnetization is 
now rotated out of its easy direction and finally aligned 



parallel to the external field if the latter exceeds a certain 
value called reorientation fieldii^i^ If the external field is 
enhanced further the magnetization does not change its 
orientation any more and the spin wave frequency v§ is 
again a linearly increasing function of the external field. 
However, before the reorientation field is reached, ^o(-Bo) 
is constant (for dipolar coupling) or decreases (for lattice 
anisotropies) as a function of the external field. There 
is an interesting feature for negative K 2 : at the reori- 
entation field the resonance frequency is zero, i.e. the 
excitation gap is closed in spite of the presence of an ex- 
ternal field and lattice anisotropies. This is not found for 
positive values of K 2 or for dipolar coupling, where the 
gap stays finite at the reorientation field (left panel). In 
such cases the probe frequency v^f may be lower than 
the resonance frequency vq for all angles A9b and no 
resonance field is detected at all. This feature is not 
reproduced by the classical model, but experimental ev- 
idence for such a scenario was indeed reported for nickel 
films with enhanced uniaxial anisotropyiii 
There are more differences between the classical and 
quantum mechanical theory for negative anisotropies and 
dipolar coupling, namely for field values that are smaller 
than the reorientation field. Consider, for example, an 
external field applied perpendicularly to the easy direc- 
tion (8b = 90°, dots). Here the resonance frequency and 
therewith the excitation gap is zero for all field values 
below the reorientation field in the classical calculation. 
The same is expected for the quantum mechanical calcu- 
lation, since in this region the magnetization has a finite 
component parallel to the film plane in the ground state 
and thus the ground state is degenerated with respect to 
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FIG. 3: The resonance fields for positive (solid lines, filled 
symbols) as well as for negative (dashed lines, open symbols) 
lattice anisotropics as a function of the polar angle of the 
external field (left panel, T = 0) and as a function of re- 
duced temperature t = T/T c (right panel, 8b = 0° - squares, 
9b — 90° - circles). The position of the resonance field of a 
paramagnet is indicated by the arrow in the right panel. 



FIG. 4: The resonance frequencies as a function of the exter- 
nal field applied in the hard direction for a coupled bilayer 
system. Left: symmetric film system, K<2, a — K^b = 5/iskG, 
right: Asymmetric film system Ki a = 5i\2i> = 5^i_BkG. For 
both panels the interlayer coupling strengths are Ji = 1 jueV 
(solid lines) and Ji — 2 /ieV (circles) . Additionally the results 
for single films, i.e. Jj = are shown in the right panel(dash- 
dotted line - layer A and dotted line - layer B). 



the azimutal angle <fi of the magnetization. Therefore a 
Goldstone mode is expected, i.e. the spin-wave energy 
for q = should be zero. This feature is not reproduced 
by our approximation, since it restricts the problem to a 
plane perpendicular to the film plane. However, despite 
these differences, the agreement between the classical and 
the microscopic theory is rather good for a large region 
of the parameter space at T = 0. 

Let us now consider the temperature dependence of the 
dispersion relation. In the following discussion the exter- 
nal field will be applied either along the easy direction or 
along the hard direction. This configurations will be re- 
ferred to shortly as "easy direction" or "hard direction" 
case. In Fig. [2] we display the dispersion relation for posi- 
tive K 2 at three different reduced temperatures t — T/T c 
ranging from t = 0.0 to t = 1.0 for the easy as well as 
for the hard direction case. Additionally the result for 
a reduced anisotropy K' 2 = §.bA7K 2 at zero temperature 
t = is shown (circles). The thin solid line denotes the 
resonance as found for a paramagnet (J a — 0). The dis- 
persion relation shows a significant temperature depen- 
dence: With increasing t the curves move towards the 
paramagnetic line. The magnitude of this shift in a cer- 
tain temperature interval is determined by the distance 
between the spin wave frequency at t = and the param- 
agnetic mode. The resonance frequencies for the easy di- 
rection case are more sensitive to temperature than those 
for the hard direction case, since the former are further 
away from the paramagnetic line. However, the param- 
agnetic solution is not reached at T c since a finite magne- 
tization is induced by the external field. Further one no- 
tices that the shape of the curves at t = and at t = 0.7 



are qualitatively the same. Temperature reduces the res- 
onance frequency for zero external field as well as the re- 
orientation field. Similar effects are found for fixed tem- 
perature but reduced anisotropy constant K 2 as can be 
seen by comparing the results for K 2 — 15 /iskG,i = 0.7 
(dashed line) and K 2 = 8.2/i B kG, t = (circles). How- 
ever, there are still differences between both cases and 
hence a reduced anisotropy parameter has clearly not 
the same effects on the dispersion relation as increasing 
temperature. The curves for temperatures in the vicinity 
or above the Curie temperature and for t = differ even 
qualitatively. Similar results are found for negative lat- 
tice anisotropies and for dipolar coupling. These findings 
bear important consequences for the usual interpretation 
of FMR data in terms of a classical model where the 
effect of temperature is taken into account via effective 
anisotropy parameters K 2 (T). Obviously, this ansatz is 
an acceptable approximation only for lower temperatures 
and not at all justified for temperatures in the vicinity 
or above T c . 

At the end of this section we want to discuss resonance 
fields (Fig. The angular as well as the temperature 
dependence is shown for positive (solid line, filled sym- 
bols) as well as for negative (dashed line, open symbols) 
lattice anisotropies. A probe frequency of 9 GHz is cho- 
sen. It is clearly seen that the resonance field has a max- 
imum at the hard direction and a minimum at the easy 
direction. The resonance fields at these extremal points 
are larger for negative than for positive anisotropies. 
If the temperature is enhanced, the resonance fields shift 
towards the paramagnetic limit i?or(param) = 2.9 kG, 
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FIG. 5: The dispersion relations for a symmetric film system 
{K 2 ( a ,b) = 5/xskG, J i — 1 fieV) for different reduced tem- 
peratures t — T/T c . The external field is applied in the hard 
direction. The paramagnetic solution is also shown (thin solid 
line) . 

which is however not reached before t « 2.5. Again the 
magnitude of these shifts in a certain temperature inter- 
val scale with the distance between the resonance field 
at zero temperature and the resonance field found for a 
paramagnet. Thus for both signs of the anisotropy the 
temperature dependence of the resonance field is more 
pronounced if the external field is applied perpendicu- 
larly to the film plane rather than parallely. The reso- 
nance fields increase with temperature in the easy direc- 
tion and decrease with temperature in the hard direction. 
All characteristic features follow from the respective dis- 
persion relations in Fig. ^and Fig. Keeping in mind 
the behavior of the single film we now want to analyze 
the coupled bilayer system to understand the intrinsic 
temperature dependence of the interlayer exchange cou- 
pling. Concerning the coupled system we will consider 
two cases: a symmetric film system where the upper and 
the lower layers are equivalent and an asymmetric film 
system, i.e. a sytem with different anisotropy strengths 
in the upper and the lower layer. 



IV. THE INTERLAYER COUPLING 

We start our discussion considering dispersion rela- 
tions. Fig. 0] shows the dispersion relation of two 
equivalent layers (left panel, K^ a = = 5/iskG) 

as well as of an asymmetric film system (right panel, 
K2a = 5/iskG, Kib = lMfikG). The interlayer cou- 
plings are 1 fieV (solid lines) and 2 /zeV (circles) , the tem- 
perature is zero. For the asymmetric system the results 
for decoupled films (Jj = 0) is also shown. The easy 
direction is perpendicular to the film and the external 
field is applied in the hard direction. For the coupled 
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FIG. 6: The same as Fig. |KJ but for an external field applied 
in the easy direction (9b q =0°). 



film two modes are found, the acoustical (lower) and the 
optical one (upper). In the case of symmetric parame- 
ters the curves are parallel and the distance between the 
modes Av scales with the interlayer coupling strength 
J i . The position of the acoustical mode is not influenced 
by the interlayer coupling strength. The same applies for 
the easy direction case. The situation is somewhat more 
complex if the layers are not equivalent as shown in the 
right panel. For Bq — the modes have lower frequen- 
cies compared to the symetric system. This is caused by 
the reduced anisotropy strenght in layer b. Furthermore 
the modes are not parallel any more and the distance be- 
tween the modes depends on both: the film asymmetry 
and the interlayer coupling. Other than for symmetric 
films the acoustical mode is now slightly influenced by 
the interlayer coupling strength Jj, It is located between 
the modes of the upper and the lower layer of the decou- 
pled system. 

Next we want to analyze the temperature dependence 
of the resonance modes at fixed anisotropy and inter- 
layer coupling parameters. Note that this is not pos- 
sible within the classical theoryji2*ii where different 
anisotropy parameters and interlayer coupling constants 
have to be assumed for every temperature. Fig. [5] and 
Fig. [S]show the dispersion relations for equivalent layers 
at three different reduced temperatures: t — 0.0, t = 0.7 
and t — 1.0 for the easy and the hard direction case as 
well as the dispersion relation of a paramagnet (thin solid 
lines) . Similar trends as for the single layer are found for 
the temperature dependence of the resonance frequen- 
cies. All modes move towards the paramagnetic solution 
with increasing temperature and the magnitudes of the 
shifts scale with the distance to the paramagnetic mode. 
For the hard direction case (Fig. [SJ) the optical mode is 
influenced only slightly in the linear regime, since the dis- 
tance to the paramagnetic mode is very small. Because 
in this regime this distance is larger for the acoustical 
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FIG. 7: The resonance fields as a function of temperature 
for the easy direction. The results of two symmetric film 
systems (Kia = K^b = 5^skG - solid lines, K^ a = K^b = 
1 /iflkG - dashed lines) as well as of an asymmetric film system 
(K2a = 5K2b ~ 5/iskG - dashed lines with circles) are shown. 
Further parameters: Ji = lfieV, v$ = 12GHz. 



mode the latter is stronger influenced by temperature. 
For zero external field the mode frequencies of the opti- 
cal and acoustical mode get smaller with increasing tem- 
perature and reach zero at T c . The distance between the 
optical and acoustical mode is clearly reduced by increas- 
ing temperature but stays finite for finite fields at the 
Curie temperature. For the easy direction case (Fig. 0) 
the t — resonance frequencies of both modes are clearly 
larger than the paramagnetic one. This causes a signifi- 
cant decrease of the resonance frequencies with increasing 
temperature. The distance between the modes Afo is re- 
duced with temperature. This reduction is comparable 
in the easy and in the hard direction case. 
Let us now consider the temperature dependence of the 
resonance fields for the coupled bilayer system. The 
probe frequency is chosen to be 12 GHz. The resonance 
field for a paramagnet is £?o r (param) = 3.87 kG. Fig. [7| 
(easy direction) and Fig. [S] (hard direction) show the res- 
onance fields as a function of temperature for ferromag- 
netic interlayer coupling (Jj = 1/ieV). Two symmetric 
films with higher (K 2 = 5 /iskG - solid lines) and smaller 
anisotropy strength (K 2 — 1 /iskG - dashed lines) as well 
as an asymmetric film system (symbols) are considered. 
In the easy direction case all resonance fields are smaller 
than the paramagnetic solution and therefore the fields 
increase with temperature. For zero temperature the res- 
onance fields are smaller for larger anisotropy strengths. 
The mode distance, however, is the same for both sym- 
metric systems (ABo r ks 1.25 kG). It is clearly seen that 
the asymmetric system is somehow a " mixture" between 
the considered symmetric systems. The acoustical as well 
as the optical resonance fields are located between the re- 
spective fields of the symmetric systems. The resulting 
distance between the optical and the acoustical field is 
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FIG. 8: The same as Fig. but for external fields applied in 
the hard direction. 

larger than for the symmetric systems. This is in agree- 
ment to the findings of Fig. Since the resonance fields 
for higher anisotropies are further away from the para- 
magnetic solution they react more sensitively to temper- 
ature than the resonance fields for the low anisotropy 
system. The distance between the modes is reduced by 
temperature for all considered systems. 
The behavior in the hard direction case is similar. The 
acoustical resonance fields are larger than the param- 
agnetic one and consequently the fields are reduced by 
temperature. The resonance fields of the optical mode 
increase with temperature for the low anisotropy and for 
the asymmetric film system since they are smaller than 
the paramagnetic field. If the resonance field is very close 
to the paramagnetic one the resonance field may behave 
non monotonically as found for the optical resonance of 
the high anisotropy system (lower solid line). As for the 
easy direction case the distance between the optical and 
acoustical resonance field A£?or is reduced by tempera- 
ture. The last point is of special interest for symmet- 
ric systems. Recall that for zero temperature as well as 
within the classical approachi2*ii the distance between 
the optical and acoustical field AB 0r is solely determined 
by the interlayer coupling strength. This feature, how- 
ever, is lost at finite temperatures. As can be seen in 
Fig- El the distance shrinks faster for the high anisotropy 
system K 2 — 5/iskG than for the system with lower 
anisotropy K 2 — I /iskG, irrespective of the same inter- 
layer coupling strength. Thus the field distance Ai?o r 
depends on both, on the interlayer coupling strength as 
well as on the anisotropies. Note that for finite tempera- 
tures this holds even for symmetric systems, i.e. systems 
with equivalent magnetic films. This behavior is not re- 
produced by the classical theor y 10 ! 11 . For the asymmet- 
ric system the temperature dependence is comparable to 
those of the high anisotropy system. Furthermore the 
temperature dependence of the distance between the op- 
tical and acoustical field ABq t is slightly enhanced for 
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FIG. 9: The distance between the optical and acoustical 
resonance field ABor as a function of the reduced temper- 
ature t — T/T c normalized to its t = value. For Ji = 
1 fieV results are shown for a symmetric system with a lower 
anisotropy (dashed line), for a symmetric system with a higher 
anisotropy (solid line) and an asymmetric film system with a 
high anisotropy layer and a low anisotropy layer (filled dots). 
For J i — 2 neV results are shown for the symmetric low 
anisotropy system (open squares). 



higher coupling strengths (open squares). 

Let us summarize our findings for the coupled bilayer 

system: 

• There are two uniform modes belonging to the 
acoustical and optical spin wave branch for the cou- 
pled bilayer system. For equivalent layers and at 
t = the mode distance Avq is a measure for the 
interlayer coupling strength. 

• Both modes and the respective resonance fields 
are shifted towards the paramagnetic solution with 
increasing temperature. The magnitude of these 
shifts is determined by the distance of the respec- 
tive mode from the paramagnetic line. The dis- 
tance between the optical and acoustical mode Ai/q 
and the distance between the resonance fields ABo r 
are reduced by temperature. 

• The magnitudes of these reductions are found to be 
larger for systems with a larger anisotropy strength. 

The temperature dependence of the resonance frequen- 
cies and fields is reflected in the temperature dependence 
of the effective interlayer coupling Ji(T) if the latter is 
obtained from experiment. It is clear that the effective 
interlayer coupling J/ (T) decreases if the temperature is 
enhanced. The magnitude of this decrease, however, may 
depend on the anisotropies of the film and is larger for 
larger anisotropies. Note that this behavior is exclusively 
due to the contribution of spin wave interactions to the 
temperature dependence of the uniform mode. 



V. CONCLUSIONS AND OUTLOOK 

Commonly there are two methods used to evaluat 
FMR experiments, namely a method based on an expan- 
sion of the free energy 11 or a method based on the equa- 
tion of motion for the classical magnetization vector. 10 
Both techniques use a classical continuum model. In this 
paper we propose an alternative method based on an ex- 
tended microscopic Heisenberg model. The theoretical 
treatment of this model was shown to give excellent re- 
sults for arbritrary directions of the external fields This 
new approach allows for an explicit consideration of tem- 
perature as far as magnetic excitations are concerned. 
We used the theory to make general predictions about 
the temperature dependence of the FMR frequencies and 
fields caused by magnetic excitations. Additionally the 
results of the new method were compared to results of 
the classical model. A single film as well as an interlayer 
coupled bilayer system were considered. We found the 
following results: 

• For zero temperature the results of the new method 
are very similar to those of the classical model for 
many parameter settings. There is however an im- 
portant difference in the hard direction: The res- 
onance frequency stays finite for all external fields 
for positive lattice anisotropies, while the classical 
model predicts a vanishing resonance frequency at 
the reorientation field. That is why the classical 
model can not explain the absence of a resonance 
field in the hard direction which was, however, ob- 
served experimentally for some systems. 1 - Using 
the microscopic analysis this observation turns out 
to be the natural consequence of the finite fre- 
quency minimum in the dispersion relation (see left 
panel of Fig. 2J. For large anisotropy parameters 
the probe frequency is lower than the minimal res- 
onance frequency and no resonance field can be de- 
tected any more. 

• For the coupled bilayer system two modes are 
found. For t = and equivalent magnetic layers the 
distance between the modes is solely determined by 
the interlayer coupling strength. For a system with 
non-equivalent layers the mode distance depends on 
the interlayer coupling as well as on the asymmetry 
of the film system. 

• The resonance frequencies and fields have a signif- 
icant temperature dependence. The dispersion re- 
lation looks qualitatively similar for temperatures 
below T c . The frequency for zero external field as 
well as the minimum of the dispersion relation are 
shifted to lower values with increasing t. This ef- 
fects are similar, but not equivalent, to a reduc- 
tion of the anisotropy parameter (see Fig. [5J. For 
temperatures in the vicinity or above T c the shape 
of the dispersion relation changes drastically com- 
pared to t = 0. 
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• With increasing temperature the resonance fre- 
quencies and fields are shifted towards the param- 
agnetic solution. This applies for the optical as well 
as for the acoustical branch and causes a decrease 
of the distance between the optical and acousti- 
cal mode At'o as well as of the distance between 
the respective resonance fields ABo r . However, the 
paramagnetic solution is not reached at T c since the 
external field induces a finite magnetization above 
T 

• For symmetric films the behavior of ABo r is note- 
worthy. At t — it is solely determined by 
the interlayer coupling strength J/. However, the 
anisotropy parameters do influence the tempera- 
ture dependence of AB 0r . The latter decreases 
faster with temperature for larger anisotropies. 
This qualitatively different behavior of the mode 
distance for zero and finite temperatures can not 
be described by the classical approach. 

These are the trends concerning the temperature depen- 
dence of the FMR signals due to the interaction of spin 
waves. 

Furthermore we expect valuable insights, if our model is 
used to analyze temperature dependent FMR measure- 



ments. Concerning the temperature dependence of the 
interlayer exchange coupling this method can be used to 
separate the different mechanisms which are known to 
be present in such systems, namely the spacer effect due 
to a smearing out of the Fermi surface within the spacer 
and the magnetic effect caused by the excitation and in- 
teraction of spin waves. The relevance of both effects is 
intensively discussed in literature but a separation and 
thus a direct comparison of both effects was not possible 
up to now. Such calculations are intended for the future. 
Additionally the procedure can be used to analyze tem- 
perature dependent FMR measurements of anisotropy 
parameters^ Again the temperature dependence due to 
spin waves can be separated from other sources of tem- 
perature dependence which opens the possibility to gain 
a deeper understanding of anisotropy effects, too. 
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